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1.  Introduction 


The  diamond  anvil  apparatus,  with  its  small  and  portable  size,  has  become  popular 
in  experimental  high-pressure  physics,  as  compared  to  the  bulky  piston-cylinder 
apparatus  (Hall  1964).  Previous  reports  (Aihaiti  and  Hemley  2008,  2010;  Huang 
2009)  have  documented  the  effort  in  using  the  diamond  anvil  technique  to 
determine  the  pressure-volume-temperature  (P-V-T)  relationship  for  the  ballistic 
gelatin,  which  is  widely  used  as  a  soft  tissue  simulant. 

The  following  steps  outline  how  to  use  the  Brillouin  scattering  technique  with  the 
diamond  anvil; 

1)  The  pressure  in  the  diamond  anvil  is  determined  by  a  standard  ruby 
fluorescence  method. 

2)  The  Brillouin  scattering  of  an  Ar-ion  laser  through  the  diamond  anvil  is 
measured  to  find  the  longitudinal  and  the  transverse  sound  velocities. 

3)  The  longitudinal  and  the  transverse  velocities  are  integrated  to  find  the 
corresponding  densities  through  the  relationship 

(P-Po)It=  ,'dP,  (1) 

where  po  and  p  are  the  density  at  pressure  Po  and  P,  respectively,  ul  and  ux 
are  the  longitudinal  and  transverse  sound  velocities,  respectively,  from  the 
Brillouin  scattering  measurement,  and  y  (=  Cp/Cv)  is  the  specific  heat  ratio. 
The  expression  joins  Ul  and  n-p  to  form  the  bulk  sound 

velocity  c,  which  relates  to  the  adiabatic  bulk  modulus  Bs  through  the 
expression  pc^  =  85.  The  specific  heat  ratio  y  turns  the  adiabatic  P-p 
relationship  data  (Eq.  1)  into  isothermal  so  that  the  P-p  relationship  data  can 
be  fitted  with  the  Vinet  or  the  Birch-Mumaghan  equations.  For  lack  of  data, 
the  ratio  y  has  been  set  to  1  following  previous  experiences  with  other 
Brillouin  measurements. 

The  earlier  primary  characterization  work  for  the  ballistic  gelatin  was  the  report  by 
Winter  and  Shifler  (1975).  The  density  is  measured  by  displacement  of  hexane.  The 
specific  heat  is  measured  by  temperature  equilibration  in  a  Dewar  flask.  In  the 
report  by  Winter  and  Shifler  (1975),  the  density  and  the  specific  heat  for  the  ballistic 
gelatin  are  reported  for  the  room  temperature  only. 

The  Brillouin  scattering  technique  provides  a  potential  of  refining  the  estimation  of 
the  bulk  modulus.  However,  it  requires  a  more  accurate  measurement  of  the  specific 
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heat  ratio.  In  this  report,  we  attempt  to  determine  the  specifie  heat  ratio  y  by  finding 
the  volume  thermal  expansion  eoeffieient  a  (from  density  variation  measurement) 
and  the  specific  heat  capacity  Cp  through  experiments  to  refine  the  P-V-T  data 
analyses.  The  P-V  data  are  fed  into  the  Vinet  and  the  Birch-Mumaghan  equations 
with  the  expectation  to  refine  the  estimation  of  the  bulk  modulus.  Furthermore,  with 
the  density,  the  sound  velocity,  the  adiabatic  bulk  modulus,  the  Griineisen 
parameter  F,  and  the  coefficient  Si  (from  literature),  a  set  of  equations  of  state 
(EoSs)  can  be  put  together  for  simulation  use. 

2.  Materials  and  Experiments 


2.1  Volume  Thermal  Expansion  Coefficient 

Ballistic  gelatin  mix  (derived  from  skin,  white  connective  tissues,  and  bones  of 
pigs;  Gelatin  Innovations  Inc.,  Shiller  Park,  IL),  rated  250  Bloom  Type  A,  was 
prepared  in  a  20%  (by  weight;  prepared  by  US  Army  Research  Laboratory  [ARE] 
Survivability/Lethality  Analysis  Directorate  [SLAB])  concentration  following  the 
NATO  standard  (Winter  and  Shifter  1975;  Nicholas  and  Welsch  2004).  The  gelatin 
blocks  for  ballistic  testing  are  prepared  and  cured  in  a  container  with  dimensions  of 
5x7x15  inches.  The  samples  are  collected  using  molds  from  the  lower  portion  of 
the  gelatin  block  (i.e.,  around  the  shot  line  during  typical  ballistic  experiments). 
The  gelatin  masses  from  the  molds  are  sliced  to  thickness  from  3  to  5  mm  (strips 
or  cubes)  to  enhance  heat  transfer  during  measurements.  Table  1  shows  the  sample 
mass  for  each  test.  The  sample  strips  or  cubes  are  loosely  filled  into  the  sample  cup 
(100  cm^),  and  then  the  sample  cup  is  gently  lowered  into  the  pycnometer  chamber 
for  measurement  (Eig.  I).  A  Micromeritics  (Norcross,  GA)  AccuPyc  II  1340 
pycnometer,  with  an  external  circulating  bath  for  temperature  control,  is  used  for 
the  measurement  of  density  as  a  function  of  temperature  from  15  to  30  °C.  Higher 
temperature  settings  are  omitted  to  avoid  phase  change  (Micromeritics  2010). 

Table  1  Sample  mass  for  each  test 


Test 

Mass 

(g) 

la 

54.0246 

lb 

50.8807 

2a 

47.8623 

2b 

35.0758 

2 


Fig.  1  Sample  cup  in  the  sample  chamber  to  be  sealed  by  the  chamber  cap  (upper  right) 


The  mass  divided  by  the  measured  volume  in  the  pycnometer  gives  the  density  for 
each  sample.  Figure  2  shows  the  calculated  density  data  along  with  their  averaged 
values.  Table  2  shows  the  measured  density  values  along  with  their  average  and 
their  standard  deviation.  The  density  at  room  temperature  reported  in  Winter  and 
Shifler  (1975)  is  approximately  1.06  g/mL. 
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Fig.  2  The  measured  density  data  along  with  their  averaged  values  (at  atmospheric 
pressure) 
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Table  2  Density,  average  density,  and  standard  deviation 


Temperatnre 

(°C) 

Density 

(g/cmh 

Density 

Density 

Density 

Average 

Density 

Standard 

Deviation 

15 

1.0655 

1.066663 

1.06599 

1.068869 

1.066756 

0.001488 

20 

1.067226 

1.068913 

1.068648 

1.072041 

1.069207 

0.002029 

25 

1.070301 

1.071428 

1.070156 

1.076338 

1.072055 

0.002911 

30 

1.077899 

1.075084 

1.074044 

1.084595 

1.077905 

0.004748 

The  volume  thermal  expansion  coefficient  a  is  evaluated  from  the  density  data 
following  the  formula; 


a  = 


(2) 


The  calculated  volume  thermal  expansion  coefficients  are  shown  in  Fig.  3.  The 
phenomenon  of  negative  thermal  expansion  has  been  observed  in  rubber  and  some 
polymers  (ref:  http://berkelevphvsicsdemos.net/node/344  [Physics  @  Berkeley, 
accessed  2015  June  30];  also  see  Miller  et  al.  2009;  Lind  2012;  Takenaka  2012; 
Fang  et  al.  2014). 
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Fig.  3  The  caicniated  voinme  thermai  expansion  eoefficients  a  from  the  density  data 
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2.2  Specific  Heat  Capacity 


The  specific  heat  capacity  is  measured  using  the  differential  scanning  calorimetry 
(DSC)  (TA  Instruments,  New  Castle,  DE)  technique  as  outlined  in  ASTM  E 
1269-05  (Robinette  2008;  Piatt  2010).  Measured  specific  heat  capacity  Cp  data  for 
samples  of  gelatin  (prepared  by  ARL/SEAD)  and  their  average  value  are  shown  in 
Eig.  4.  Table  3  shows  the  measured  specific  heat  capacity  values  along  with  their 
average  and  their  standard  deviation.  The  specific  heat  measured  at  room 
temperature  reported  in  (Winter  1975)  is  approximately  1.13  cal/g/°C 
(=  4.73  J/g/K). 
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Fig.  4  The  measured  specific  heat  capacity  Cp  data  and  the  average  value  (line  marked  “rl” 
came  from  Rohinette  2008;  lines  marked  “pi”,  “p2”,  “p3”  came  from  Piatt  2010) 


Table  3  Specific  heat  capacity,  average  heat  capacity,  and  standard  deviation 


Temperature 

(°C) 

Cp 

|J/(gK)| 

Cp 

Cp 

Cp 

Average 

Cp 

Standard 

Deviation 

20 

4.787 

5.099 

4.736 

5.411 

5.008 

0.313 

25 

5.364 

5.073 

4.706 

5.372 

5.129 

0.314 

30 

6.109 

5.069 

5.018 

5.445 

5.410 

0.503 
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3.  Data  Analysis 


3.1  The  Cp/Cv  ratio  (=  y) 

From  the  measured  specific  heat  capacity  at  constant  pressure  (Cp)  data  and  the 
volume  thermal  expansion  coefficient  (a)  data,  the  specific  heat  capacity  at  constant 
volume  Cv  can  be  calculated  using  the  relationship 

(Cp  -  C,)  =  ^  Bt  ,  (3) 

where  Cp  is  the  specific  heat  capacity  at  constant  pressure,  Cv  is  the  specific  heat 
capacity  at  constant  volume,  a  is  the  volume  thermal  expansion  coefficient,  T  is  the 
temperature,  p  is  the  density,  and  Bt  is  the  isothermal  bulk  modulus.  From  the 
calculated  Cv  the  ratio  of  Cp/Cv  can  be  calculated.  Figure  5  shows  the  Cp/Cv  ratio 
(sometimes  called  the  adiabatic  index)  in  the  temperature  range  tested  using  the 
averaged  values  from  measured  data  and  a  typical  isothermal  bulk  modulus  of 
2.6  GPa  (estimated  from  the  measured  sound  speed  of  0.156  cm/ps  [Winter  and 
Shifler  1975]). 
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Fig.  5  The  Cp/Cv  ratio  in  the  temperature  range  tested 
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3.2  The  Gruneisen  Parameter 


The  Mie-Griineisen  EoS,  as  widely  applied  in  simulation  codes,  uses  a  Gruneisen 
parameter  to  relate  to  the  thermal  pressure 


r  =  V 


a  Bj  _  a  Bj 

CvP  Cpp  ^  ’ 


(4) 


where  E  is  the  Gruneisen  parameter,  a  is  the  volume  thermal  expansion  coefficient, 
and  Y  is  the  Cp/Cv  ratio.  Eigure  6  shows  the  Gruneisen  parameter  versus 
temperature. 


< 

-0.05 

9  C 

<  <■ 

3  i' 

■«<  r 

4  < 

4  ( 

t  « 

«  ( 

)  0 

4  <■ 

9  C 

4  r 

3  r 

•4 

Q)  -0.15 

0^ 

re  -0  2 

ra 

^  -0.25 

c 

0) 

► 

❖ 

c 

o 

ej 

< 

► 

temperature  (“C) 

Fig.  6  The  Gruneisen  parameter  vs.  temperature 

Negative  Gruneisen  parameter  has  been  reported  for  amorphous  ice  (Andersson  and 
Inaba  2005),  amorphous  silicon  (Eabian  and  Allen  1997),  graphene  (Zakharchenko 
et  al.  2009;  Bera  2011),  germanium  (Harris  and  Avrami  1972),  and  carbon  dioxide 
in  overdriven  explosion  products  (Medvedev  2014).  The  Gruneisen  parameter  and 
the  thermal  expansion  coefficient  have  been  found  to  display  a  sign  change  close 
to  a  quantum-critical  point  (Garst  and  Rosch  2005).  Analysis  of  thermal  expansion 
of  a  classical  chain  with  pair  interactions  performing  longitudinal  and  transverse 
vibrations  shows  that  the  thermal  expansion  coefficient  and  the  Gruneisen 
parameter  can  be  nonlinear  and  negative  (Kuzkin  and  Krivtsov  2015).  The  negative 
sign  of  the  Gruneisen  parameter  is  probably  associated  with  its  transverse  mode  but 
not  its  longitudinal  mode  (Andersson  and  Inaba  2005;  Kuzkin  and  Krivtsov  2015). 
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4.  Fitting  Equations 


4.1  The  Vinet  EoS 

The  calculated  Cp/Cv  ratio  (=  y)  is  used  in  Eq.  1  to  find  the  P-V  data,  which  are 
fitted  into  the  Vinet  EoS, 


where  Bxo  is  the  isothermal  bulk  modulus  at  zero  pressure  Po,  B^o  is  its  pressure 
derivative  at  zero  pressure  Po,  and  Vo  is  the  initial  specific  volume  at  zero  pressure- 
to  find  the  isothermal  bulk  modulus  Bxo  and  its  pressure  derivative  B^o  at  zero 
pressure. 

Since  the  measured  longitudinal  and  transverse  velocities  show  a  discontinuity  at 
approximately  2  GPa  (Pig.  7)  and  a  phase  change  has  been  identified  under  the 
microscope  near  that  pressure,  only  the  data  below  2  GPa  are  used  in  the  data  fitting 
procedure. 
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Fig.  7  Measured  longitudinal  and  transverse  velocities  from  the  Brillouin  scattering  tests 
(from  Aihaiti  and  Hemley  2010) 
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Figure  8  shows  the  fitted  eurves  for  various  derivatives  B^o  of  the  isothermal  bulk 
modulus  and  the  data  from  the  Brillouin  tests  (line  marked  DAC). 


Fig.  8  Fitted  curves  to  the  Vinet  equation  of  state  along  with  the  data  from  the  Brillouin 
tests 

If  the  fitting  algorithm  is  set  free  to  find  both  the  isothermal  bulk  modulus  Bto  and 
its  derivative  B^o,  the  values  found  are  B^o  =  16.91  and  Bto  =  3.801,  whieh  is  higher 
than  expeeted.  If  the  derivative  B^o  is  preseribed,  various  values  of  the  isothermal 
bulk  modulus  B  are  found  and  listed  in  Table  4; 


Table  4  Fitted  values  for  the  Vinet  equation 


Derivative  B^o 

Isothermal  Bulk  Modulus  Bto 

14 

4.345 

0.9943 

16 

3.963 

0.9957 

16.91 

3.801 

0.9958 

20 

3.291 

0.9942 

22 

2.996 

0.9916 

24.82 

2.6 

0.9859 

26 

2.479 

0.9832 
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The  sensitivity  of  the  fitting  procedure  is  lower  than  expected. 


4.2  The  Birch-Murnaghan  EoS 


Another  widely  used  equation,  the  Birch-Mumaghan  EoS,  is  also  exercised  to  fit 
the  data: 


(P-Po)It 


(6) 


where  Bto  is  the  isothermal  bulk  modulus,  B^o  is  its  pressure  derivative  at  zero 
pressure,  and  Vo  is  the  initial  specific  volume  at  zero  pressure. 

Figure  9  shows  the  fitted  curves  for  various  derivatives  B^o  of  the  isothermal  bulk 
modulus  and  the  data  from  the  Brillouin  tests  (the  line  marked  DAC). 


density  (g/cm^) 


Fig.  9  Fitted  curves  to  the  Birch-Murnaghan  EoS  along  with  the  data  from  the  Brillouin 
tests 

If  the  fitting  algorithm  is  set  free  to  find  both  the  isothermal  bulk  modulus  Bto  and 
its  derivative  B^o,  the  values  found  are  B^o  =  27.28  and  Bto  =  3.248,  which  is 
higher  than  expected.  If  the  derivative  B^o  is  prescribed,  various  values  of  the 
isothermal  bulk  modulus  Bto  are  found  and  listed  in  Table  5. 
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Table  5  Fitted  values  for  the  Birch-Murnaghan  equation 


Derivative  Bjg 

Isothermal  Bulk  Modulus  Bxo 

24 

3.508 

0.9958 

27.28 

3.248 

0.996 

36 

2.713 

0.9953 

38.21 

2.6 

0.995 

40 

2.522 

0.9947 

The  sensitivity  of  the  fitting  proeess  is  lower  than  expeeted. 

4.3  Comparison  with  the  Shock  Hugoniot  Compression  Data 

Shoek  Hugoniot  eompression  data  using  fiat  plate  impact  experiments  on  gelatin 
(Nagayama  et  al.  2006;  Shepherd  et  al.  2009;  Appleby- Thomas  et  al.  2011)  can  be 
summarized  in  these  Us  -  Up  relationships: 

Us  =  1.57  +  1.77  Up  ,  (7) 

Us  =  1.45  +  1.99  Up  .  (8) 

Both  of  these  equations  are  converted  into  the  Steinberg’s  equation  (without  the 
energy  term) 

[l-(Si-l)p]2  ’ 

where  p,  =  (—  —  1 )  is  the  compression,  and  Si  comes  from  the  Us  -  Up  relationship. 

These  shock  Hugoniot  data  are  plotted  along  with  the  adiabatic  data  from  the 
Brillouin  scattering  measurements  as  shown  in  Fig.  10.  The  omitted  energy  term 
involves  integration  over  the  time  of  deviatoric  stresses  and  pressure,  which  is 
difficult  for  desktop  arithmetic.  However,  the  plotted  line  is  very  close  to  the  P-V 
data  found  in  the  paper  by  Appleby- Thomas  et  al.  (2011),  which  shows  that  the 
omitted  energy  term  is  very  small  in  the  plotted  pressure  range. 
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Fig.  10  The  Brillouin  scattering  adiabatic  data  and  the  shock  Hngoniot  data 

The  discrepancy  between  the  Brillouin  data  and  the  shock  Hugoniot  data  hints 
toward  probable  inaccuracy  associated  with  the  Brillouin  scattering  data  for  this 
soft  material.  While  further  investigations  are  still  needed  to  find  the  exact  causes, 
the  inaccuracies  may  come  from  the  conditions  that  the  soft  material  may  no  longer 
be  homogeneous  or  isotropic  within  the  diamond  anvil  during  the  tests  and  that 
there  is  a  difference  in  the  strain  rates  (quasi-static  for  the  diamond  anvil,  while  a 
high  strain  rate  for  the  shock  Hugoniot);  and  the  measurement  may  also  be  under 
the  influence  of  some  traces  of  deviatoric  components  emerging  during  the  tests. 
Furthermore,  the  bulk  moduli  calculated  from  the  measured  longitudinal  and 
transverse  velocities  (ul  and  u-r)  (Fig.  7)  also  appear  to  be  higher  than  expected. 

5.  Summary  and  Conclusions 

The  ballistic  gelatin  samples  are  measured  to  find  their  density  and  the  specific  heat 
capacity  as  a  function  of  temperature.  The  volume  thermal  expansion  coefficient, 
the  specific  heat  capacity  ratio  (from  tests  in  the  range  of  20-30  “C),  and  the 
Griineisen  parameter  as  a  function  of  temperature  are  calculated  from  the  measured 
data.  These  values  are  used  to  improve  the  P-V  relationship  data. 

These  P-V  data  are  then  fitted  to  the  Vinet  and  the  Birch-Mumaghan  EoS. 
However,  the  fitted  curves  show  low  sensitivity  of  the  equations  to  the  derivatives 
of  the  bulk  modulus.  Moreover,  the  P-V  data  from  the  Brillouin  scattering 
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measurement  show  substantial  differences  compared  to  the  published  shock 
Hugoniot  data.  The  discrepancies  with  the  shock  Hugoniot  data  hinted  toward 
inaccuracies  in  the  current  measured  Brillouin  scattering  data  (the  longitudinal  and 
the  transverse  velocities)  for  this  soft  material. 

However,  with  the  density,  the  sound  velocity,  the  adiabatic  bulk  modulus,  the 
Griineisen  parameter  T,  and  the  coefficient  Si  (from  the  literature),  a  set  of  EoS  can 
be  put  together  for  simulation  use. 
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Appendix.  Synthetic  Polymer  Gels 
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A,1  Synthetic  Polymer  Gels 


The  US  Army  Research  Laboratory  (ARL)  is  actively  researching  the  thermo¬ 
plastic  elastomer  gels  (TPEGs),  which  is  a  class  of  nonaqueous,  synthetic  gels  for 
ballistic  testing.  TPEGs  are  typically  composed  of  a  triblock  copolymer  with 
glassy  or  crystalline  end  blocks  separated  by  a  rubbery  midblock  and  swollen  with 
a  midblock-selective  solvent.  The  mechanical  properties  of  the  material  can  be 
altered  through  the  incorporation  of  a  solvent  that  will  selectively  swell  the  rubbery 
midblock. 

The  poly(styrene-b-ethylene-co-butylene-b-styrene)  (SEBS)  triblock  copolymer 
and  the  mineral  oil  as  solvent  mixed  at  the  ratio  32.5%  (volume)  SEBS  and  67.5% 
mineral  oil  produces  a  synthetic  gel  with  compatible  ballistic  performance.  It  is 
tested  for  its  volume  thermal  expansion  coefficient  for  comparison.  Table  A-1 
shows  the  sample  mass  for  each  test.  Figure  A-1  shows  the  measured  density  as  a 
function  of  temperature.  Table  A-2  shows  the  measured  density  values  with  their 
average  and  their  standard  deviation,  along  with  the  density  of  the  ballistic  gelatin 
(average,  cf  Fig.  2  in  this  report)  for  reference.  The  density  in  the  temperature  range 
is  lower  than  that  of  water. 

Table  A-1  Sample  mass  for  each  test 


Test 

Mass 

(g) 

ml 

31.9689 

m2 

30.0406 

m3 

29.6690 

^  Mrozek  RA,  Leighliter  B,  Gold  CS,  Beringer  IR,  Yu  JH,  VanLandingham  MR,  Moy  P,  Foster  MH, 
Lenhart  JL.  The  relationship  between  mechanical  properties  and  ballistic  penetration  depth  in  a  viscoelastic 
gel.  J  Mech  Behavior  Biomed  Materials.  2015;44: 109-120. 

2 

Foster  M,  Moy  P,  Lenhart  J,  Mrozek  R,  Weerasooriya  T.  Punch  Response  of  Gels  at  Different  Loading 
Rates.  Aberdeen  Proving  Ground  (MD):  US  Army  Research  Laboratory  (US);  2014  Mar.  Report  No.:  ARL- 
TR-6882.  Also  available  at:  http://www.arl.amiv.mil/arlreports/2014/ARL-TR-6882.pdf 
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Fig.  A-1  Measured  density  data  with  the  density  for  the  hallistic  gelatin 


Table  A-2  Density,  average  density,  and  standard  deviation 


Temperature 

(°C) 

Density 

(g/cmh 

Density 

Density 

Average 

Density 

Standard 

Deviation 

15 

0.8742 

0.873 

0.8738 

0.873667 

0.000611 

20 

0.8718 

0.8707 

0.8716 

0.871367 

0.000586 

25 

0.8693 

0.869 

0.8696 

0.8693 

0.0003 

30 

0.8671 

0.8669 

0.8671 

0.867033 

0.000115 

35 

0.8639 

0.8662 

0.865 

0.865033 

0.00115 

Figure  A-2  shows  the  volume  thermal  expansion  coefficient  calculated  from  the 
measured  density  variations  with  the  thermal  expansion  coefficient  of  the  ballistic 
gelatin  (average,  cf.  Fig.  3)  for  reference.  The  average  value  can  be  taken  to  be  a 
constant  at  approximately  0.0005 
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Fig.  A-2  The  calculated  volume  thermal  expausiou  coefficieuts  from  the  deusity  data 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ARL 

US  Army  Research  Eaboratory 

DSC 

differential  scanning  calorimetry 

EoS 

equation  of  state 

P-V-T 

pressure-volume-temperature 

SEBS 

styrene-b-ethylene-co-butylene-b-styrene 

SLAB 

Survivability/Eethality  Analysis  Directorate 

TPEG 

thermo-plastic  elastomer  gel 
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